Autonomy and Singularity in Dynamic Fracture 
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The recently developed weakly nonlinear theory of dynamic fracture predicts 1/r corrections to 
the standard asymptotic linear elastic l/\/r displacement-gradients, where r is measured from the 
tip of a tensile crack. We show that the 1/r singularity does not automatically conform with the 
notion of autonomy (autonomy means that any crack tip nonlinear solution is uniquely determined 
by the surrounding linear elastic 1/s/r fields) and that it does not automatically satisfy the resultant 
Newton's equation in the crack parallel direction. We show that these two properties are interrelated 
and that by requiring that the resultant Newton's equation is satisfied, autonomy of the 1/r singular 
solution is retained. We further show that the resultant linear momentum carried by the 1/r singular 
fields vanishes identically. Our results, which reveal the physical and mathematical nature of the 
new solution, are in favorable agreement with recent near tip measurements. 
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A weakly nonlinear theory of dynamic fracture, which 
extends the standard theory of fracture [l], [4], was re- 
cently developed [3, |4[ . This theory was shown to be in 
excellent agreement with groundbreaking experimental 
measurements of the deformation near the tip of rapid 
cracks [3|-l6J. Furthermore, it may be relevant to under- 
standing currently poorly understood crack tip instabili- 
ties [7[ . In this Rapid Communication we derive a series 
of theoretical results that further elucidate the physical 
and mathematical nature of the theory. 

The standard approach to dynamic fracture - linear 
elastic fracture mechanics (LEFM) [lj, |2J - assumes that 
linear elasticity dominates the deformation fields outside 
a small region near the tip of a crack. Its major prediction 
is that crack tips concentrate large deformation-gradients 
and stresses which are characterized by a universal l/-\/F 
singular behavior. Many results in this theoretical frame- 
work are derived from the latter property [jj, |2j ■ 

The basic physical idea underlying the weakly nonlin- 
ear theory of dynamic fracture is that linear elasticity 
breaks down when elastic nonlinearities intervene near 
the tip of a crack. This is a physically intuitive idea since 
atomes/molecules are expected to sample reversible (i.e. 
elastic) anharmonic parts of the interaction potential be- 
fore their separation is large enough to induce irreversible 
deformation (e.g. damage, plasticity and eventually frac- 
ture). 

To mathematically formulate this idea, we consider the 
following expansion of the displacement field u [3, |j] 

u(r, 9; v) ~ eu (1) (r, 0; v) + eV 2) (r, 9; v) + C(e 3 ) 

= w (1) (r, 9; v) + u i2) (r, 9; v) + C(e 3 ) , (1) 

where e quantifies the magnitude of the displacement- 
gradients [3] , {r, 9) is a polar coordinates system located 
at a crack's tip and moving with it at a speed v in the 
9 = direction. The first order term in e corresponds to 
linear elasticity, which is actually only a first term in a 
more general expansion, while the second order term cor- 
responds to the leading nonlinearity that intervenes when 
the deformation is large enough. As higher order nonlin- 
earities are neglected in Eq. ([!]), the theory based on it 



is termed "weakly nonlinear theory of dynamic fracture" 

nan. 

The expansion in Eq. (Q]) can be substituted in a gen- 
eral elastic strain energy functional U(F), where Fy = 
Sij + djUi, from which the first Piola-Kirchhoff stress ten- 
sor s can be derived as 



dU 

dF 



(2) 



s quantifies forces in the deformed configuration per unit 
areas in the reference configuration [9|. Then, the mo- 
mentum balance equations and the crack faces traction- 
free boundary conditions are obtained order by order in 
e dHll]. M (1) ( r >#5«) in Eq. ([I]) satisfies the first order 
problem, which is a standard LEFM one [1, 2]. The near 
crack tip (asymptotic) fields for steady state propagation 
under Mode I (opening) symmetry are given by 
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u£\r,6;v) = ^Ln x (0:ri. 



tW(r,9-v) = 
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n y (9;v). 



(3) 



Here Kj is the Mode I "stress intensity factor" which 
cannot be determined by the asymptotic analysis, but 
rather by the global crack problem. $1(0; v) is a known 
universal function jl|-l3] and /j, is the shear modulus, x 
corresponds to the propagation direction, 9 = 0, and y 
corresponds to the direction in which the tensile loadings 
are applied, # = ±7r/2. The displacement fields in Eq. ([3]) 
give rise to the famous l/\/r displacement-gradients and 
stress singularity |l|, |2j ■ 

u^ 2 \r, 9; v) in Eq. ((T|) satisfies the second order prob- 
lem, which was explicitly derived in [3, |j|. It has the 



following form 
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A, / v 2 sin 2 9 
A log r + — log I 1 ^ — 



i?a s log H — — log I I 



v 2 sin 2 I 
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Aa^d-B^ + Tj,^;^ 



tan 9 d , s = ay, s tan 0, a 2 s = 1 - v 2 /c^ s 



(4) 



where Cd ]S are the dilatational and shear wave speeds, 
respectively. Y(0;u) is given in the form 



T x (6; v) = 2_c n {v) cos(n0), 



T y (0; v)=J~)i n (v) sin(n0), 



where the coefficients {c n (v), d n (v)} can be easily ob- 
tained by solving a set of linear algebraic equations [3|, |j] . 
The coefficients A and B are related by the traction-free 
boundary conditions on the crack faces [3j, |j] through 
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2/iBa s — (A + 2fi)dgT y (ir; v) — ^jlk{v) 
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Here 



k(v) = — 



(A + 2[i)a 2 



16a 2 v 4 (A + /i) 



Hci[Aa s a d -(l + a 2 s ) 2 f 



(6) 



(7) 



and A is the second Lame coefficient [l|, [2j . The displace- 
ment fields in Eq. @ contain log (r) terms and give rise 
to 1/r singular displacement-gradients. Both of these 
features were directly verified by near tip measurements 

ia. 

The solution in ((U)-© satisfies the second order mo- 
mentum balance equations and the traction-free bound- 
ary conditions on the crack faces. Therefore, one may 
reach the conclusion that the second order solution con- 
tains a parameter B that is not uniquely determined by 
the stress intensity factor Kj. If true, this result has 
profound theoretical implications as it suggests that the 
concept of the autonomy of the near crack tip nonlinear 
region [J Q is not always valid. The basic idea behind 
this concept is that the mechanical state within the near 
tip nonlinear zone, which is surrounded by the LEFM 
fields of Eqs. (|3]), is uniquely determined by the value of 
Ki and is otherwise independent of the applied loadings 
and the geometric configuration in a given problem. This 
implies, for example, that systems with the same Ki, but 
with different applied loadings and geometric configura- 
tions, will be in the same mechanical state within the 
near tip nonlinear zone. Autonomy is a central concept 
in fracture mechanics |l|, |2j . 

In |3j values of B were directly extracted from the 
experimental data, without addressing the question of 
whether they can be in fact theoretically determined. In 
[J] it was shown that in the quasi-static limit, v— > 0, B 



can be theoretically determined by Ki and hence the au- 
tonomy of the near tip nonlinear region is retained. This 
may appear as a puzzling result, because Eq. (j4|), with 
Eq. (J6]), satisfies the asymptotic second order boundary- 
value problem for all B's - what is then the missing physi- 
cal ingredient that is not contained within the asymptotic 
boundary- value problem? 

To answer this question, which was only partially ad- 
dressed in [4], consider the net force per unit sample 
thickness / acting on a line of radius r encircling a crack's 
tip 



f, 



SijTijrd9 



(8) 



where n is an outward unit normal on the circle. The 
(5) Mode I symmetry immediately implies f y = 0. Moreover, 



as no force is acting on the crack tip in the crack propa- 
gation direction (x), and focusing first on the quasi-static 
limit (v — > 0) in which material inertia does not play a 
role, we must have f x — 0. In the framework of quasi- 
static LEFM, one can show that a 1/r contribution to 
s results in f x =^ 10]. Therefore, f x — is satisfied 
within quasi-static LEFM if and only if the 1/r singular- 
ity is discarded altogether, i.e. its prefactor is set equal 
zero [lOj. The 1/r singularity generates an unbalanced 
(spurious) force in the crack parallel direction (where no 
boundary conditions are imposed), even though it satis- 
fies the asymptotic boundary-value problem. Hence, in 
the framework of LEFM this singularity is unphysical. 

This conclusion is in sharp contrast to the correspond- 
ing situation in the quasi-static weakly nonlinear theory 
[J]. In this case, a 1/r singular contribution to s also 
does not automatically lead to f x = 0, but the latter can 
be recovered without discarding the whole solution by 
properly choosing B as a function of Kj in Eqs. ((4]) and 
([6|), retaining autonomy |4|. Therefore, the special prop- 
erty of the 1/r singularity discussed above is precisely 
the missing physical ingredient, which is not contained 
within the the asymptotic boundary-value problem, that 
ensures that autonomy is not violated. We thus see that 
in contrast to LEFM, the 1/r singularity in quasi-static 
weakly nonlinear fracture mechanics is a physically sound 
solution that does not violate any physical principle and 
conforms with the concept of autonomy. 

The discussion above was restricted to the quasi-static 
limit. The generalization to the fully dynamic case, v > 0, 
is somewhat more subtle because material inertia can 
play a role. To understand this, we write down the resul- 
tant Newton's equation for the material enclosed within 
a circle of radius r around the tip. It is a balance be- 
tween the force / in Eq. ([8]) and the time rate of change 
of linear momentum p (both per unit sample thickness) 

Sijrijrd9 — v 2 p I r'dr' I d xx Uid6=pi, (9) 

-7T JO J-TT 

where the steady state relation dt = —vd x was used and 



r' is a dummy integration variable. We expect Eq. (|9]) 
to provide the necessary condition for determining B in 
the dynamic case, though it is clear that /^0 is possible 
if p 7^ 0, without violating any physical law. 

In order to demonstrate the latter possibility, consider 
the asymptotic first order solution given in Eq. ([3|) , which 
can be used to derive the standard first order (linear elas- 
tic) stress tensor s^ 1]. Using u^ and s^ in Eq. (|9|), 
/W and p^ can be calculated. Recall that Mode I sym- 
metry implies that f y — p y — 0, so all the discussion to 
follow focuses on the the crack parallel direction x. 
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FIG. 1. (Color online) r" 1 / 2 /! 1 '/^/ vs. v/c s . An identi- 
cal curve is obtained when r~ Y ' 2 p x jKi is plotted vs. v/c s , 
which shows that f x = p x ■ 



In Fig. Uwe plot r" 1 ! 2 f^ jKi vs. v/c s . The fig- 
ure shows that f x ^ for v > (except for an iso- 



lated point). Moreover, a direct calculation shows that 
fx —Px ■ Therefore, the standard LEFM l/y/r singu- 
larity automatically satisfies Eq. (j9)), though fx y= 0. 
The situation is different when the 1/r singularity of the 
weakly nonlinear theory is considered. 

Equation Q was shown above to be satisfied automat- 
ically to first order in e, when the fields in Eq. ^ arc 
used. Consider now Eq. ([9]) to second order in e. The 
second order stress tensor s^ has the following scaling 
property 



^-du^du^+du^ 



(10) 



where the tensorial notation was omitted for simplicity. 

(2) 

This implies that fx is a, constant independent of the 
radius r. On the other hand, we have 



d x 



■(2) 



,( 2 ) 



(11) 



which implies that p x depends on r (in fact it diverge 
logarithmically). Therefore, the only way in which Eq. 
([9]) can be satisfied to second order is by having px = 0, 
which implies that 



/. 



(2) = 



,(2) 



s x jn jr d6 = 



(12) 



Therefore, the 1/r singular fields carry no net linear mo- 
mentum and f x =0 is the condition that determines 5, 
precisely as in the quasi-static limit. 

(2) 

The latter prediction, i.e. that p x =0, can be directly 
checked using the explicit second order solution in Eq. 
(|3]). In order to calculate p x , we evaluate the different 
contributions to d r 



,(2) 



d X x^x(0;v) = ^-n c n (v)sin9 : 



i cos(n8) sin 6 + 2 cos 6 sin(n6) 
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Alogr + -logl 1- — -J 



v 2 sin 2 I 



v 2 sin 2 i 



Ba s 
Ba s log r + — — ^ log I 1 - 



2Ac 2 d [v 2 + (2c 2 , - v 2 ) cos(26>)] 2Bc 2 s ^/l - v 2 /c 2 [v 2 + (2c 2 s - v 2 ) cos(20)] 



[2c 2 - v 2 + v 2 cos(29)} r 2 



[2c 2 - v 2 + v 2 cos(26>)] r 2 



(13) 
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Hence, calculating analytically the angular integrals over we obtain 
the above expressions, which sum up to J_ d xx u x dd, 

I 



d x 



A\ogr+ 2-log{ I 



Asin(26) 



4cos(20) 



v 2 sin 2 I 



Ba s 
Ba s log r H — - log 1 



Ba s sin(26) 



l + a 2 + l 4 cos(20) 



v 2 sin 2 i 



T x (8;v) 



d9 



(14) 



+ y^ nc n (v) 



sin[(n - 2)0] - 2sin(n0) + sin[(n + 2)6"] 

At 2 



= 



which explicitly verifies that the 1/r singular fields carry 
no net linear momentum, p x =0. 
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FIG. 2. (Color online) d y u y (r,0) of the weakly nonlinear the- 
ory (solid line), cf. Eq. fl]), together with the experimental 
data (circles) of 0, 0| for v = 0.53c s and the material used in 
|5J]. The theoretical curve is obtained without any fit, where 
only the stress intensity factor reported in Fig. 1(b) of [3| was 
used as an input, see text for more details. 



We reiterate that surprisingly, even at finite crack ve- 
locities (v > 0) inertia does not play a role in determin- 
ing B and in retaining autonomy, and the condition to 
be satisfied remains as in the quasi-static limit, i.e. Eq. 
(fT2| . As was stressed several times above, Eq. (fT2]) is 
not satisfied for any B, but rather determines it. There- 
fore, in order to calculate B, we substitute s^ (which 
is obtained from expanding Eq. (J2J in orders of e, cf. 
[3, |4| ) in the integrand of (fT2j) and look for the value of 
B that makes the integral vanish. This cannot be done 
analytically, but is easily achieved numerically. 

To test the theory, we compare its predictions to 
the direct near-tip deformation measurements of [5j. 
For that aim, we focus on v = 0.53c s and use 
{c n (0.53c s ),d n (0.53c s )}, \ — 2^jl, \i = 32. 5kPa as reported 
in [3J for the material used in [5| . Using these numbers in 
Eqs. ©-© t° obtain s^ and then calculating numeri- 
cally the integral in Eq. (112ft . we obtain B ~ 18.5. This 
value cannot be directly compared to the value obtained 
by a fitting procedure in Fig. 1(b) of [3j because in the 



latter case the asymptotic LEFM fields of Eqs. {3} had to 
be supplemented with a subleading term (corresponding 
to the "T-stress" [3J, |5j) and hence an additional param- 
eter was involved. Instead, we focus on a smaller region 
near the tip, where the subleading term that we do not 
consider here is less significant, and use the stress inten- 
sity factor of Fig. 1(b) in [|, Ki = 1250Pa v / ^- Recall 
that, as required by autonomy, the stress intensity factor 
is the only parameter that is needed as an input to the 
asymptotic near tip theory. In Fig. [2] we plot d y u y (r,0) 
of the weakly nonlinear theory, cf. Eq. ([T]), together with 
the experimental data of |3j, [5J . The agreement between 
the theory and the experimental data is remarkable, sup- 
porting the theoretical results derived above. 

In summary, in this Rapid Communication we theo- 
retically explored some properties of the 1/r singularity 
in the framework of the recently developed weakly non- 
linear theory of dynamic fracture. It was shown that the 
theory is consistent with the notion of the autonomy of 
the nonlinear near-tip region for any crack tip velocity 
v, extending the quasi-static results of [J]. In addition, 
it was shown that no net linear momentum is carried 
by the 1/r singular fields. As only Mode I symmetry 
was considered here, a direction for future investigation 
is the development of a weakly nonlinear theory for more 
general fracture conditions. 



[1 
[2 
[3 
[4. 
[5 

[e; 

[7] 
[8] 

[9 
[10 



L. B. Freund, Dynamic Fracture Mechanics, (Cambridge 

University Press, Cambridge, 1998). 

K. B. Broberg, Cracks and Fracture, (Academic Press, 

1999). 

E. Bouchbinder, A. Livne and J. Fineberg, Phys. Rev. 

Lett. 101, 264302 (2008). 

E. Bouchbibder, A. Livne and J. Fineberg, J. Mech. Phys. 

Solids 57, 1568 (2009). 

A. Livne, E. Bouchbinder and J. Fineberg, Phys. Rev. 

Lett. 101, 264301 (2008). 

A. Livne, E. Bouchbinder, I. Svetlizky and J. Fineberg, 

Science 327, 1359 (2010). 

E. Bouchbinder, Phys. Rev. Lett. 103, 164301 (2009). 

E. Bouchbinder and T. S. Lo, Phys. Rev. E 78, 056105 

(2008). 

G. A. Holzapfel, Nonlinear Solid Mechanics, (Wiley, 

Chichester, 2000). 

J. R. Rice, J. Mech. Phys. Solids 22, 17 (1974). 



